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Abstract. We study the resonant tunneling properties of an electron through a few types of binary periodic 
and aperiodic multibarrier systems. Within the framework of the effective-mass approximation, we calculate 
the transmission coefficients to investigate the dependence of the transmission resonances on the system 
parameters such as the kind of aperiodicity, the generation number, and the widths of the wells and 
barriers. Similarities and differences of the resonances between the binary periodic and aperiodic systems 
are discussed in detail. Transmission resonances in aperiodic systems are found to be characterized by 
complex resonance splitting and a variety of peak-to- valley ratios which are not exhibited in the periodic 
system. For some energy ranges, transmission resonances in aperiodic systems are also found to resemble 
those in the periodic system, despite the existence of aperiodicity. 

PACS. 73.40.Gk Tunneling - 71.55.Jv Disordered structures; amorphous and glassy solids - 73.20.Dx 
Electron states in low-dimensional structures 

1 Introduction beam epitaxy and metal-organic chemical vapor deposi- 
tion, there has been a lot of work on the problem of the 

Tunneling of an electron through a potential barrier is electronic resonant tunneling in semiconductor superlat- 

one of the fundamental phenomena in quantum mechanics tices [2-15]. Among interesting features emerged from the 

and plays a key role in the physics of electronic and opto- studies, the most well recognized features are the reso- 

electronic devices [1]. Stimulated by the advancement in nance splitting effects and the energy band effects; for 

modern material fabrication technology such as molecular- a finite superlattice composed of N identical potential 
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(N — l)-fold resonance splitting, and the split resonant en- 
ergies approach the band structure for large N [1, 12-14]. 
Very recently, Guo et al. [15] studied the resonance split- 
ting effects in superlattices which are periodically juxta- 
posed with two different building barriers to demonstrate 
that the resonance splitting is determined not only by the 
structure but also by the parameters of building barriers. 

In a different context, there has been much interest 
in the electronic properties of deterministic aperiodic sys- 
tems [16-19] which are known to have more complex ge- 
ometrical structures than the periodic system. Studies on 
these systems have revealed a variety of exotic electronic 
properties such as the singular continuity of the energy 
spectrum, the self-similarity of the electronic wave func- 
tions, the power-law behavior of the resistance, and so 
forth. However, most of the work has focused on the elec- 
tronic properties in the infinite limit of the system size 
[16-20], and the study on the transport properties of the 
finite-size systems, particularly from tunneling point of 
view, has received less attention. Recently, Singh et al. 
[21] studied the electronic transport properties of the Fi- 
bonacci and Thue-Morse (TM) superlattices to compare 
the results with those of the periodic system. Besides, Liu 
et al. [22] calculated the electronic transmission spectra 
of the Cantor fractal multibarrier systems to show that 
the tunneling spectrum is more complex than that of the 
periodic system. 

In this paper, we study the aperiodicity-induced ef- 
fects on the transmission resonances in a few types of 
binary aperiodic multibarrier systems. To do this, tak- 
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Fig. 1. Schematic configuration of a part of the multibarrier 
systems used in the calculation. 

ing into account four kinds of binary multibarrier systems 
whose geometrical structures arc determined by Eq. (11), 
we calculate the transmission coefficients of an electron 
through these multibarrier systems. Dependence of the 
transmission resonances on the system parameters such 
as the kind of apcriodicity, the generation number, and 
the widths of the wells and barriers is investigated. From 
this, we first illustrate the characteristics of the trans- 
mission resonances exhibited in the binary periodic (BP) 
system, and then make a comparison of the resonances 
between the BP and aperiodic systems; similarities and 
differences between them are presented in detail. In doing 
this, a comparison of the resonance splitting exhibited in 
the common-well (CW) structure with those exhibited in 
the common-barrier (CB) structure of the systems is also 
presented. 

2 Method 

We shall now derive an expression for the transmission 
coefficient of a multibarrier system using the transfer ma- 
trix formalism. To do this, we consider an electron with a 
longitudinal energy E incident from left of the system. As- 
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suming that the phonon scattering can be neglected and jth and (j + l)th wells as 
no bias is applied across the system, the effective-mass ap- 
proximation leads to the continuous Schrodingcr equation 
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(6) 



2m* dx 2 



V(x) 



$ j (x)=E* j (x) 



(1) 



where Tj is the unimodular transfer matrix, and oij and 



where x is the longitudinal direction of the system, V(x) 
the minimum energy of the conduction band, and to* = 



l w(b) 



the effective mass of the electron in the well (bar- 



rier) of the jth cell. Figure 1 shows the schematic config- 
uration of a part of the system. For convenience of calcu- 
lation, we set the conduction band minimum to be zero 
and the potential barrier to be rectangular, i.e., 



V(x) 



for Xj < x < yj 



(2) 



V for yj < x < Xj+i 

where Xj (yj) is the starting position of the jth well (bar- 
rier). Then, the wave function associated with the electron 
in the jth cell can be written as 



^•(a;) = A 3 e ik( - X - X ^ + B,e- ik( - x ~ x ^ 



(3) 



in the well and 



i(x) = Cje- K{x - yi) + Dje K(x - %) 



(4) 



in the barrier. Here, 



2m* E 



l2m* b (V -E) 



(5) 



/3j are given by 
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Here, Wj(= yj — Xj) and bj(— Xj+i — yj) are the widths of 
the jth well and barrier. Multiplying Tj successively, we 
can write the relation of ^4's and -B's between the first and 
the (N + l)th region as 



A N+1 \ / A 1 

B )" MN \b 

where M/v is the total transfer matrix given by 



(8) 



Mpf — TjvTjv_i • • • T 2 Ti = 



qn b 



N 



b* N a* N 



(9) 



Since there will be a reflected wave in the first region but 
only a transmitted wave in the (N + l)th region (i.e., 
_Bjv + i = 0), we can write the transmission coefficient as 

f 



T 



(10) 



l Nl 



are the wave numbers in the wells and barriers, respec- 
tively. 

By applying the Bastard's matching conditions of the 
wave function and its derivative at discontinuity of V(x), 
we can write the relation of the coefficients between the 



Generally, it requires extensive matrix manipulation to 
calculate T. However, for the multibarrier systems con- 
sidered in this paper, T can be easily calculated in terms 
of the deterministic substitution rules given below. 

We now introduce four kinds of deterministic sequences 
- the BP, the TM [16,17], the period-doubling (PD) [18], 
and the copper-mean (CM) [19] sequences which are gen- 
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erated by the substitution rules with a\ = ola&b + P\Pb and b\ = olbPa + ol* a Pb for the 

BP sequence, 

BP: S l+1 =Sf, S 1= AB , 

_ _ a;+i = aiai + 6;&f , k +1 = aibi + a*bi 

TM : S l+1 = SiS tl (S ,S ) = (B,A) <_ _ _ (14) 

( n ) a ;+1 = aiai + bib[ , b l+1 = aj&j + affy 



PD : S l+1 = StSl,, (So^S,) = (A,AB) 
CM : S l+1 = SiSf_ lt (S_i,S ) - (B,A), 



with ao = a a, &o = Pa, o-o = Qg, and 60 = Pb for the TM 
sequence, and 



ai+i = Maf-! + l&Z-i| 2 ) + fo ; *^-i( a '-i + a *-i) 

(15) 

b l+1 = 6i(af_ a + |&;-i| 2 ) + a,*&;_i(a ; _i + af_ x ) 
with ao = o-a, bo = /3a, cl\ = Q-aQ-b + P\Pb, and b± = 
&bPa + ®*aPb for the PD sequence. Recursion relations 
for the CM sequence are exactly the same as Eq. (15) with 
a_i = a B , b-i = Pb, ao = &A, and b = Pa- 

3 Numerical results and discussion 



where I is the generation number (i.e., N — 2 l for the 

first three sequences and TV = [2 l+2 — (— 1)']/3 for the 

last sequence) and Si is the complement of Si which is 

obtained by exchanging the letters A and B. Here, A and 

B represent two kinds of unit cells of the multibarrier 

system. The unit cell A (B) consists of a well with the 

width wa {wb) and a barrier with the width bA (&s) and 

the height V. We refer the case of wa — wb with bA i=- bs 

as the CW model and the case of bA — bs with wa ^ Wb 

as the CB model, respectively As a sample material for calculation, we choose the /xc- 

By means of Eq. (11), we can write the recursion rela- Si:H/a-Si:H supcrlattice [22], where ^c-Si:H acts as the 

tions of the total transfer matrices between different gen- well and a-Si:H the barrier with V = 0.4 eV. The effective 

orations as mass in the wells and barriers is taken to be m* w = ml = 

0.3m e [23], where m e is the free electron mass. In cal- 
culating transmission coefficients of an electron through 
(12) multibarrier systems, we treat two cases separately; the 
one is to set the widths of the wells equal while arrange 
the widths of the barriers according to the given substi- 

TT . . . . . ^ ,,„. .... , tution rule (the CW model), and the other is to set the 

Using the relations m Eq. (12), we can easily derive the 

. widths of the barriers equal while arrange the widths of 

recursion relations of as and s between different gener- 

. ,, the wells according to the given substitution rule (the CB 

ations as follows: 

model) . Some examples of the results are plotted in Fig- 

&;+i = a i + b*bi ure s 2 and 3. In plotting, the mesh of E is taken to be 

(13) 

bi+i = h(ai + a ; *) AE = 1.0 x 10~ 5 eV, and the system parameters are set 



BP : Mj+i = Mf, Mi = T B T A 

TM : Mj+i = M1M1, (M , M ) - (T B ,T A ) 
PD : M i+1 = Mf_ x M u (Mo, Mr) - (T A ,T B T A ) 
CM : M +1 = Mf_ x M u (M_i,M ) - {T B ,T A ). 
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to be w A = w B = 20 A and b A = 2b B = 8 A in the CW ° 

model, and b A = b B = 7 A and wa = 2w B — 40 A in 
the CB model, respectively. Figure 2 shows the results ob- 
tained from the CW model for the energy range near the 
lowest domain of resonances, and Figure 3 shows the re- 
sults obtained from the CB model for the energy range 
near the first two lowest domains of resonances. Here, the 
'domain of resonances' means the energy range that con- 
tains resonant peaks in the finite-size system and would 
approach the allowed energy band in the infinite limit of 

the system size. 

0.06 0.14 0.22 0.06 0.14 
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Fig. 2. log 10 T(E) versus E near the lowest domain of reso- 
nances for the BP [(a) and (e)], TM [(b) and (f)], PD [(c) and 
(g)], and CM [(d) and (h)] systems with wa — wb — 20 A and 
bA — 2bs = 8 A (the CW model). The generation number is 
I = 2 for (a)-(d) and I = 5 for (e)-(h). 



We first discuss the features of the transmission res- 
onances exhibited in the BP system with I = 2. In this 
case, three resonant peaks in each domain of resonances 
are expected due to overlap of quasi-bound states in the 
three well regions, and the result obtained from the CW 
model [Figure 2a] agrees well with the expectation. An 

interesting feature to be noted in Figure 2a is that the the second lowest domain are suppressed. The second is 
first and the third peaks of the three resonant peaks are that the resonance widths exhibited in the CB model are 
complete (i.e., T = 1) while the second peak is incomplete much narrower and sharper than those in the CW model. 
(i.e., T < 1). We will see that the two complete peaks We now discuss the features of the transmission res- 
locate in the middle of the subdomains while the incom- onanccs in the BP system with large I. As I increases, 
plete peak disappears for large I [see Figure 2e] . As for the successive resonant splitting effects and the energy band 
result obtained from the CB model [Figure 3a], there are effects on the transmission properties are expected, and 
two distinctive features from the result obtained from the we confirm them. Figure 2e shows the transmission coeffi- 
CW model. The one is that there occurs suppression of cients obtained from the CW model with I — 5. Here two 
resonant peaks. We can see in Figure 3a that there exist distinctive features from the case of the periodic system 
a single complete peak in the first lowest domain and two [12-14] are emerged. The first is that the main domain 
complete peaks in the second lowest domain, which implies splits into two subdomains, the centers of which corre- 
that two peaks in the first lowest domain and one peak in spond to the first and the third resonant peaks in Fig- 
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Fig. 3. log 10 T(E) versus E near the first two domains of res- 
onances for the BP [(a) and (e)], TM [(b) and (f)], PD [(c) 
and (g)], and CM [(d) and (h)] systems with wa = 2wb — 40 
A and bA = &s = 7 A (the CB model). The generation 
number is I — 2 for (a)-(d) and I = 5 for (e)-(h). 



different behavior; no splitting into subdomains occurs, 
which indicates that the binary periodicity of the system 
does not affect on the splitting of this domain. 

We also study the features of the transmission reso- 
nances with varying the widths of the wells and barriers, 
and observe that the number of resonance domains in- 
creases with increasing the width of the well, which can be 
easily understood by noting that the energies of the quasi- 
bound states in a well are approximately in proportional 
to the inverse square of the width of the well [24] . We also 
observe that the widths of resonance domains decrease 
and the peak-to-valley ratios increase with increasing the 
width of the barrier. 



Having seen the features of the transmission resonances 
in the BP system, we now discuss the features of the trans- 
mission resonances exhibited in the aperiodic TM, PD, 
CM systems with I = 2. In this case, the three, three, and 
ure 2a. We argue that the occurrence of this kind of reso- four resonant peaks are expected to exist in each domain 
nance splitting attributes to the binary periodicity of the of resonances due to overlap of quasi-bound states in the 
system. The second is that the total number of resonant three, three, and four wells of the TM, PD, CM systems, 
peaks is not exactly the same as the number of the wells. The results obtained from the CW model for these sys- 
In Figure 2e, each subdomain contains 15 peaks and the terns are plotted in Figures 2b — 2d, where it can be seen 
total number of peaks are 30 while there are 31 wells in that the number of resonant peaks fits with the expecta- 
the system with I — 5. The result obtained from the CB tion. However, the results obtained from the CB model 
model with I = 5 is plotted in Figure 3c. Similar features [Figures 3b — 3d] do not always fit with the expectation: 
appeared in the CW model are also exhibited in the sec- The transmission coefficients for the second lowest domain 
ond lowest main domain; the domain splits into two sub- agree with the expectation, while the number of resonant 
domains and the total number of resonant peaks is less peaks in the first lowest domain is less than the expected 
than the number of the wells. However, the transmission number due to suppression of the peaks. The number of 
coefficients in the lowest main domain display somewhat suppressed peaks is two, one, and one in the TM, PD, CM 
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system, respectively. As for the effects of aperiodicity, we eventually compose of Cantor-like fractal structure in the 
would like to mention two points. The first is that the low- infinite limit of the system size. The second is that, despite 
est resonant peaks shift towards the lower energy region, the existence of aperiodicity, the transmission resonances 
compared with that of the BP system. The second is that for some energy ranges resemble those exhibited in the BP 
there coexist the complete and incomplete resonant peaks, system. For example, from 0.101 eV to 0.128 eV and from 
For the peaks exhibited in the CW model, all the peaks 0.184 eV to 0.188 eV in Figure 2f, the transmission cocffi- 
in the TM system and the second peak in the CM system cients exhibit the 'resonance plateaus' as in the case of the 
are complete, while all the peaks in the PD system and all BP system. We also find that, for higher energy ranges, the 
the peaks expect the second peak in the CM system are effect of aperiodicity weakens such that the number and 
incomplete. In the meanwhile, for the peaks exhibited in the widths of resonant plateaus increase, which resembles 
the CB model, all the peaks are incomplete, which implies the feature of the periodic system. The results obtained 
that it is more difficult for an electron to tunnel through from the CB model of the TM, PD, and CM systems [Fig- 
the CB structure than the CW structure. Here it should ures 3f — 3h] reveal similar features to those exhibited in 
be noted that the incomplete resonant peak appeared in the CW model; complex resonance splitting effects and 
the BP system [Figure 2a] and the incomplete peaks ap- the existence of the resonance plateaus are clearly seen, 
peared in the aperiodic systems show different behavior An example of the resonance plateaus can be seen in Fig- 
as I increases; the former fails to locate in the middle of ure 3g with the energy range from 0.178 eV to 0.197 eV. 
the subbands and disappears, while the latter survives to 
locate in the middle of the subdomains 4 Summary 

We now discuss the features of the transmission reso- In summary, we studied the effects of aperiodicity on the 

nances in the aperiodic systems with large I. Figures 2f — transmission resonances of an electron through a few types 

2h show the results obtained from the CW model of the of binary aperiodic multibarrier systems with finite sys- 

TM, PD, and CM systems with I = 5, which reveal two tern size. Dependence of transmission resonances on the 

distinctive features to be mentioned. The first is that, even system parameters was investigated, from which the sim- 

though the aperiodicity of each system is binary as in the ilarities and differences of the resonances between the bi- 

BP system, the resonance splitting pattern is much more nary periodic and aperiodic systems were presented. In 

complex than that of the BP system such that a variety doing this, a comparison of resonance splitting exhibited 

of peak-to-valley ratios appear. Wc argue that the level in the common-well structure with those exhibited in the 

of hierarchy in the resonance splitting will become deeper common-barrier structure of the multibarrier systems was 

with the increase of I and the split resonant peaks will also made in detail. It was found that complex resonance 
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splitting and a variety of peak-to-valley ratios which are 
not exhibited in the periodic system are emerged as a re- 
sult of introducing aperiodicity. It was also found that 
the transmission resonances for some energy ranges in the 
aperiodic systems resemble those in the periodic system, 
despite the existence of aperiodicity. We hope that ex- 
otic resonant tunneling properties of the binary aperiodic 
multibarrier systems such as the complex resonance split- 
ting effects, deep levels of hierarchy in the peak-to-valley 
ratios, and the existence of tunneling plateaus can be ap- 
plied in designing a new type of electronic device. 
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